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Abstract. A G-strand is a map g : (f, s) 6 1 x 1 -> g(t, s) £ G into a Lie group G 
that follows from Hamilton's principle for a certain class of G-invariant Lagrangians. Some 
■ G-strands on finite-dimensional groups satisfy 1 + 1 space-time evolutionary equations that 

admit soliton solutions as completely integrable Hamiltonian systems. For example, the 
S'0(3)-strand equations may be regarded physically as integrable dynamics for solitons on 
a continuous spin chain. 

We study collisions of solutions with singular support (e.g. pcakons) on Diff (R)-strands. 
Here Diff(M) is the group of diffcomorphisms of the real line K, for which the group product 
is composition of smooth invertiblc functions. In the simplest case of peakon-antipcakon 
collisions, the solution reduces to solving Laplace's equation and is written in terms of 
Q\ ■ harmonic functions in the parameters (t,s). We also consider a complexified system of G- 

\& | strand equations for G = Diff(R) corresponding to a harmonic map g(z, z) : C —> Diff(R) 

and find an explicit expression for its peakon-antipeakon solutions. Instead of the Laplacian, 
the peakon-antipeakon solutions in the complexified case involve the wave operator. 



Key words: Hamilton's principle, continuum spin chains, Euler-Poincare equations, Sobolev 
norms, singular momentum maps, diffcomorphisms, harmonic maps 
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1 Euler-Poincare equations for a G-strand 

The Euler-Poincare (EP) theory of G-strands is an extension of the classical chiral models 
[131 [Li] to include norms that are not bi- invariant. The classical chiral models are nonlinear 
relativistically invariant Lagrangian field theories on group manifolds. As such, they are fun- 
damental in theoretical physics. The vast literature of results for these models is fascinating. 
For example, it is well known that these models are integrable in 1 + 1 dimensions and possess 
soliton solutions. See [21 [231 ED E21 El [22], [191 GO] and references therein for discussions of the 
many aspects of integrability of the chiral models, including the famous dressing method for 
explicitly deriving the soliton solutions of these models, which is given in [32]. The solitons 
for the 0(3) chiral model are particularly familiar [30] . because this model allows an integrable 
reduction to the one-component sine-Gordon equation (e.g. see [22]). Many generalizations of 
these models have been introduced. For example, generalized chiral models with metrics that 
are not ad- invariant on the Lie algebra are considered in |25j . Other generalizations of chiral 
models for non-semisimple groups are studied in [llj . An integrable chiral model in 2+1 dimen- 
sions was proposed in [28] and the possibility for fermionic interpretation of the current variables 
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was explored in [29]. Finally, their formulation as Euler-Poincare equations for G-strands with 
Lagrangians given as norms on tangent spaces of finite-dimensional Lie groups was given in 

nana. 

The class of problems treated here generalizes the classical chiral model, which appear as a 
special case of the G-strands when their structure group is SO(3) and the norm defining their 
Lagrangian is bi-invariant, as discussed in p3]. G-strands are also found here to be associated 
with a certain generalization of harmonic maps into Lie groups. Harmonic maps into Lie groups 
were studied mathematically, and their relation to the chiral model via complexification of its 
independent coordinates was discussed in [27]. Following this same idea, we shall complexify 
the independent coordinates (s,t) of the G-strand on the Lie group of diffeomorphisms Diff(R) 
and thereby derive a class of equations which may be regarded as partial differential equations 
for generalized harmonic maps of C into Diff(K). 

Main content of the paper. 

Section [T] derives the class of G-strand equations studied here in the context of Euler-Poincare 
theory. 

Section [2] discusses singular (peakon) solutions and their collision dynamics for the G-strand 
equations with G = Diff(R). 

Section [3] introduces the complexified Diff (R)-strand equations and determines their solutions 
corresponding to peakon collision dynamics. 

Left G-Invariant Lagrangian. We begin with the following ingredients of EP theory. For 
more details and discussion, see [T3| [T7], 

• Let G be a Lie group. A map g(t, s) : R x R — > G has two types of tangent vectors, u g 
and v g , corresponding to its two types of derivatives, with respect to s and t. 

• Assume that the function L(g,u g ,v g ) : TG x TG — > R is right G-invariant. 

• Right G-invariance of L permits us to define I : g x g — > R by 

l{u g g~ l iVgg" 1 ) = L(g,Ug,v g ). 

Conversely, this relation defines for any I : g x g — > R a right G-invariant function L : 
TG x TG — > R. 

• For a map g(t, s) : R x R — > G define 

u(t,s) := g t g~ l (t,s) and 7(t, s) := 5 S 5 _1 (*, *)■ 

Lemma 1. The right-trivialized tangent vectors v(t,s) and "f(t,s) at the identity of G satisfy 

7t ~ = - ad^7 . (1) 

Proof. The proof is standard and follows from equality of cross derivatives gt s = g s t, cf. [13^117]. 
As a consequence, equation (pQ) is often called a zero- curvature relation. ■ 

Theorem 1 (Euler-Poincare theorem). 

With the preceding notation, the following four statements are equivalent ]1S\ 
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i Hamilton's variational principle on TG x TG 

5 I 2 L(g(t,s),g t (t,s),g s (t,s))dsdt = (2) 

holds, for variations 5g(t,s) of g(t,s) vanishing at the endpoints in t and s. 

ii The function g(t, s) satisfies Euler-Lagrange equations for L on G, given by 

dL d dL d dL 
dg dt dgt ds dg s 

iii The constrained variational principl^ 

5 l(v(t,s),-f(t,s))dsdt = (4) 
Jti 

holds on g x g, using variations of v and 7 of the forms 

8v = Ct — ad^C and £7 = ( s — ad 7 C , (5) 
where £(t, s) £ g vanishes at the endpoints. 
iv The Euler Poincare equations hold on g* x g* 
d 51 _ 51 d 51 _ 51 n 

^ + ad ^ + ^ + ad ^=°- (6) 

Evolutionary G-strand. We now define the fundamental quantity of interest in the remainder 
of the paper. 

Definition 1. A G-strand is a map R x R — > G from R x K into a Lie group G, whose 
dynamics in (t,s) G R x R may be obtained from Hamilton's principle for a G-invariant reduced 
Lagrangian I : g x g — > R, where g is the Lie algebra of the group G. The G-strand system of 
evolutionary partial differential equations for a left G-invariant reduced Lagrangian consists of 
the zero- curvature equation ([TJ and the Euler- Poincare (EP) variational equations ©. 

Subclasses of the G-strand maps contain the principal chiral models of field theory in theoret- 
ical physics, reviewed, e.g., in [2D]. An interpretation of the G-strand equations as the dynamics 
of a continuous spin chain is given in [13] . This is the source of the term, 'strand'. The cor- 
responding theory of molecular strands is discussed [3]. Recently, a covariant field theory of 
G-strands in higher dimensions for G = Diff(R) has also been developed [5]. 



Overview and organization of the paper. In section[2j we shall develop G-strand dynamics 
for diffeomorphisms on the real line and show that the Diff(R)-strand dynamics admit singular 
solutions associated with a pair of momentum maps. In section [3] we introduce the complexified 
Diff (R)-strand equations and solve their partial differential equations corresponding to two- 
peakon collision dynamics. Finally, Section [4J summarizes the present results and provides some 
outlook for future research. 

As with the basic Euler-Poincare equations [21], this is not strictly a variational principle in the same sense 
as the standard Hamilton's principle. It is more like the Lagrange d'Alembert principle, because we impose the 
stated constraints on the variations allowed [IT] . The prool of equations ^ follows the same pattern as the proof 
of the zero-curvature relation fl]). 
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2 The G-strand PDE for G = Diff(R) 

This section studies the G-strand system arises when we choose G = Diff (R) and the Lagrangian 
involves the H 1 Sobolev norm. This case is reminiscent of fluid dynamics and may be written 
naturally in terms of right-invariant tangent vectors (t, s, x) and j(t, s, x) defined by 

dt9 = vog and d s g = 7 o g , (7) 

where the symbol o denotes composition of functions. In this right-invariant case, the G-strand 
system of partial differential equations (PDE) with reduced Lagrangian £(1/, 7) takes the following 
form, which generalizes to R rf in any number of spatial dimensions, 

d_^ + d_^_ ad *^_ ad *^ 

dtbv ds bj v bv 7 ^7 ' ^ 

d"/ dv 

The distinction between between the maps (^,7) : R X R — > q x g and their pointwise values 
(V(i, s),7(t, s)) € g x g will always be clear in context, so that no confusion will arise. Likewise, 
for the variational derivatives bijbv and b£/bj. 

Equations ([5]) form a subset of the equations studied in [f)J [T2J HS1 ES] for complex fluids. They 
also form a subset of the equations for molecular strands studied in [3]. The latter comparison 
further justifies the name G- strands for the systems being studied here. These equations may 
be derived from Hamilton's principle for an affine Lie group action, under which the auxiliary 
equation for 7 may be interpreted as an advection law. This interpretation is discussed further 
in [3J|S]. The covariant form of these equations has been developed in [5]. 

2.1 The G-strand Hamiltonian structure 

Upon setting m = bijbv and n = bl/bj, the right-invariant G-strand equations in (|S]) for maps 
R x R — t- G = Diff (R) in one spatial dimension may be expressed as a system of two 1+2 PDEs 
in (t, s, x), 

m t + n s = - ad*m - ad*n = -(um) x - mv x - {^n) x - wy x , 
7t - v s = - ad 7 i/ = -wy x + jv x . 

The corresponding Hamiltonian structure for these Diff (R)-strand equations is obtained by 
Legendre transforming to 

h(m,i) = (m, v) -£(u, 7) . 

One may then write the m-7 equations (J9|) in Lie-Poisson Hamiltonian form as 



d 


m 




di 


.7. 





— ad^m d s + ad* 
d s — ad 7 



bh/bm = v 
bh/bj = — n 



(10) 



where the symbol □ denotes {sA^m)v := ad*m for v € Q and m € Q* . This is the Lie-Poisson 
bracket dual to the action of the semidirect-product Lie algebra 

q = X(R)(S)A 1 (Dens)(R) C{d B ) , 

in which 36(R) is the space of vector fields and A 1 (Dens) (R) is the space of 1-form densities on 
the real line R, plus a generalized 2-cocycle C(d s ). In the 7 — v entry of the Hamiltonian matrix 
in (|10p . one recognizes the vector-field covariant derivative in s, and finds its negative adjoint 
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operator in the v — 7 entry. The Lie-Poisson bracket in (jlOp first arose in the Hamiltonian 
formulation of chromohydro dynamics, i.e., the dynamics of a Yang- Mills fluid plasma [5]. 
It then re- appeared in the theory of complex fluids. For discussions of how such Lie-Poisson 
Hamiltonian structures arise in complex fluids with finite-dimensional order-parameter groups 
(broken symmetries) and full details of their Lie algebraic properties, see [6J [T2"l [T5j [26] and 
references therein. 



Relation to the Camassa-Holm equation pQ. An interesting subcase of the system of 
Diff (R)-strand equations (|8]) arises when one chooses the Lagrangian ^(^,7) in (|8|) to depend 
only on u, in particular, as its H 1 norm on the real line, 

t(v,l) = l\H 2 Hl, (11) 

with vanishing boundary conditions, as \x\ — > 00. In this case, m = v — u xx , and the corre- 
sponding restriction of the equations in ([8]) provides an extension of the completely integrable 
Camassa-Holm (CH) equation [I], 

81 

m% = — ad*m = — (um) x — mu x with m = —— = v — v xx . (12) 

Specifically, these modified G-strand equations reduce in the absence of 7-dependence precisely 
to the CH equation, which admits singular solutions known as peakons in the form [1] 

m(t,x) = Y / M a(t)8(x-Q a (t)), (13) 

a 

where we sum in a € Z over the integers, or over any subset of the integers. The peakon solution 
(|13p of the CH equation has been understood as a singular momentum map obtained from the 
left action of Diff (K) on embeddings of points on the real line [TB]. The generalization to higher 
dimensions is also possible. See [16] for further discussion. 



2.2 Peakon solutions of the Diff (R)-strand equations in ([9]) 

With the following choice of Lagrangian using the H 1 norm, 

^7) = ^ll*-^il7||^, (14) 

the Diff (R)-strand equations ([9]) admit peakon solutions in both momenta m and n, with con- 
tinuous velocities v and 7. We state this result in the following theorem. 

Theorem 2. The Diff (WL)-strand equations (0) admit singular solutions expressible as linear 
superpositions summed over a G Z \1J$ 

m. 



(s, t,x) = J2 M a (s, t)5(x - Q a {s, t)) , 

a 

(s, t,x) = J2 N a (s, t)8(x - Q a (s, t)) , 

a 

u(s, t,x) = K *m = ^2 M a( s , t)K(x, Q a ) , 

a 

7(s, t, x) = K * n = N a (s, t)K(x, Q a ) , 



n{ 

(15) 
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that are peakons in the case that K(x,y) = l x V V These singular solutions follow from 
Hamilton's principle 5S = for the constrained action S = J L(v, 7, Q) alt given by 

S= I £(u, ~f) + Yl M «( s ' *) (W(*> t) ~ "(Q a i *, tj) + N a(s, t) {d s Q a (s, t) - 7 (Q a , s, t))ds dt . 

•* a a 

and they are peakons in (s,t) for the Lagrangian ^(^,7) given in equation 

The solution parameters {Q a (s, t), M a (s, t), N a (s, t)} with a € Z that specify the singular 
solutions (|15p are determined by the following set of evolutionary PDEs in s and t, in which we 
denote K ab := K(Q a , Q b ) with integer summation indices a, b,c,e € Z: 

d t Q a (s, t) = u(Q a , s,t) = J2 M b(s, t)K ab , 

b 

d s Q a (s, t) = 7 (Q a , s,t)=J2 N b (s, t)K ab , 

b 

\ dK ac (16) 

d t M a (s, t) = - d s N a - Y,(MaM c + N a N c )-— (no sum on a), 

c 

BK ec 

d t N a (s,t) = d s M a + ^(iV 6 M c - M 6 iV c )^— (K eb - K cb ){K~ l ) ae . 

b,c,e 



The last pair of equations in (116ft may be solved as a system for the momenta, or Lagrange 
multipliers (M a ,N a ), then used in the previous pair to update the support set of positions 
Q a (t,s). Given Q a (t,s) for a G Z, one constructs (m,n(t, s, x)) along the solution paths x = 
(J a (i, s) from the first row of (|15|) and then obtains (1/, 7(4, s, x)) for iGK from the second row. 
Alternatively, knowing the position Q a (s,t), a € Z, for all s at a given time i, also determines 
N a upon inverting the matrix K ab in the second equation in (|16|) . 

Single-peakon solution. The single-peakon solution of the Diff (M)-strand equations (J9]) is 
straightforward to obtain from (|16p . Namely, writing the first two equations for a single peakon 
and substituting into the third equation shows that Q 1 (s,t) satisfies the Laplace equation, 

(% + %)Q 1 (s,t)=0. (17) 

Thus, any harmonic function h(s,t) provides a solution Q 1 = h(s,t). One then finds from the 
first two equations in (fT6|) that 

M\s, t) = -^ht(s, t) N\s, t) = -^h s (s, t), (18) 

where Kq = K(0). The solution for the single-peakon parameters Q 1 ,M 1 and iV 1 depends only 
on the choice of the harmonic function h(s,t), which in turn depends on the (s,t) boundary 
conditions. In particular, the single-peakon does not depend on the shape of the Green's function 
K(x). 

2.3 Two-peakon interactions on G-strands 

Let the Lagrangian be expressed via an operator norm 

= ^\\u\\ 2 Q = 1 fuQudx, (19) 
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where the operator Q is symmetric, positive-definite and spatially translation-invariant, with 
vanishing boundary conditions, as \x\ — > oo. For example, the H 1 norm for the CH peakons in 
the Lagrangian (|14p arises when Q is chosen to be the Helmholtz operator, so that Qu = u — u xx . 

Consider peakons at positions Q 1 (s, t) and Q 2 (s, t) on the real line. Since the Green's function 
K for the operator Q depends only on the difference of the peakon positions, one naturally defines 
the new variable X(s,t) = Q 1 — Q 2 . Then the first two equations in (]16p imply 

d t X = (M 1 -M 2 )(K -K(X)), 
d s X = (N 1 -N 2 )(K Q -K(X)), 

where Kq = K(0) is the value of the Green's function K{X) for the operator Q when the peakon 
positions coincide, so that X = 0. 

The second pair of equations in (|16p may then be written as 

d t h'h = -dsNx - (MiM 2 + iViiV 2 )A"(X) , 
8 t M 2 = -d s N 2 + (M X M 2 + N X N 2 )K'(X) , 

d t N x = d s M x + (AqM 2 - M X N 2 ) K J~ K K\X) , ( 21 ) 

K -K 



8 t N 2 = 8 S M 2 + (JViM 2 - M\N 2 ) 



K + K 

Asymptotically, when the peakons are far apart, the system (|2ip simplifies, since k° +k — > 1 
and K'(X) ->■ as \X\ -> oo. 

The system ()2ip has two immediate conservation laws obtained from their sums and differ- 
ences, 

d t {M x +M 2 ) = - dsim + N 2 ) , 

d t {N 1 -N 2 ) = d s (M 1 -M 2 ). 

These may be resolved by setting 

Ml -M 2 = J^-, N!-N 2 



Kq-K' Kq-K 1 (23) 

M l + M 2 = d s <P, Aq + N 2 = - d t <j) , 

and introducing two potential functions, X and (f>, for which equality of cross derivatives will 
now produce the system of equations (j20j) and (|2T|) . Namely, from (|23l) . one may solve for 

Ma = , &tX - + -d a <f>, N t = , 9sX , - ^ t 0, 
2(A -A) 2 2(K -K) 2 W ' 

M 2 = ; r H — a s ra , A 2 = ; 0+0 ■ 

2(K -K) 2 ^' 2(K -K) 2 V 

Now substituting into the first and third equations of the system (|2ip yields the determining 
equations for the potentials X and <^>, 

2(Ao + A-) (25) 



where one denotes A = d 2 + d 2 and V = (d s ,dt). 
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A simplification. A simplification arises if (f> = 0, in which case the collision is perfectly 
antisymmetric, as seen from equation (|23|) . This is the peakon-antipeakon collision, for which 
the second equation in (f25j) reduces to 



A - Y + ^¥) |vx|2 = - <26> 

This equation can be easily rearranged to produce Laplace's equation: 

r x 

AF(X)=0, where F(X) = (K - K(Y))~ 1/2 dY . (27) 

Jx 

Thus, the dynamics of the relative spacing X(s, t) = Q 1 — Q 2 = 2Q in the peakon-antipeakon 
collision may be obtained by elementary means. In particular, for the H 1 peakon case when 
K(Y) = ie"l y l, we have 

F(X) = [ - 1 dY = cosh" 1 (e~ m/2 ) = h(s,t), (28) 
Jxo y/l - e -\ y \ V 1 

where h(s, t) is any harmonic function on the (s, t) plane. Consequently, the motion of X follows 
from 

X(s, t)=± In ( sech 2 (/i(s, t))) , (29) 
for any given choice of the harmonic function h(s,t). 

3 Complexified G-strand Diff(R) equations 

As mentioned in the Introduction, the class of problems being treated here generalizes the 
classical chiral model, which is a special case of the G-strands when their structure group is 
SO (3), as discussed in [H]. Harmonic maps are also related to the chiral model as discussed in 
[27j . upon complexifying the independent coordinates of the chiral model. Pursuing this idea, 
we shall complexify the independent coordinates (s, t) of the G-strand equations ([9]) on the Lie 
group of diffeomorphisms Diff (R) and thereby derive a class of equations which may be regarded 
as partial differential equations for harmonic maps of C into Diff (R). This is potentially a rich 
area for further mathematical study. 

3.1 Complexifying (s,t) for a real Lagrangian 

Upon complexifying (s, t) € M 2 — > (z, z) G C where z denotes the complex conjugate of z and 
setting 7 = v the Euler-Poincare G-strand equations in ([9]) become 

d_M + d_^_ ad *^_ ad *^ 

dz 5u dz 5v h ' 8v "Sv' /oq\ 
dv dv 

a = ad ^ • 

oz oz 

Here the Lagrangian £ is taken to be real, and we will assume that it is defined in terms of a 
suitably chosen norm on complex functions 



1 
2 



i.^ u ) = oIMIq = o / vQvdx, (31) 
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where Q is a real, symmetric, positive-definite operator, invariant under spatial translations in 
and with vanishing boundary conditions as \x\ — > oo. 

Since the Lagrangian I is chosen to be a suitable norm on complex functions defined by 
the spatial differential operator Q, we may regard the first equation in the complexified G- 
strand system (I30p as describing geodesic motion according to the map (z,z) £ C — > Diff(R) 
with respect to the Q-norm. The second equation (|30p may be regarded as a compatibility (or 
zero-curvature) condition in the complex parameters z and z. 

Upon setting m = 5£/8u, m = 51/59, for the real Lagrangian I, equations (f30l) may be 
rewritten as 

m z + rhz = - ad*m - ad^m = -(um) x - mv x - (um) x - rhu x , 
- Vz = ~ adpu = -vv x + vv x , 

where the independent coordinate x € R is on the real line, although coordinates (z, z) € C are 
complex, as are solutions v, and m = Qv, by equations ([30]) for the real Lagrangian (f3T|) . 



Remark 1. 



• Equations (|32|) are invariant under two involutions, P and C, where 
P : (x, m) — > (—x, —m) and C : Complex conjugation. 



• For real variables m = m, v = v =: u and real evolution parameter z = z =: t, equations 
(132|) reduce to the Camassa-Holm equation (fT2|) . 



3.2 Peakon solutions of complexified Diff (M)-strand equations in ( 13 2 j) 

The complexified Diff (R)-strand equations (j32[) admit peakon solutions in complex momenta m, 
with continuous complex velocity v. We state this result in the following theorem. 

Theorem 3. The complexified Diff(R) -strand equations $3ty) admit singular solutions expressed 
as linear superpositions summed over a € Z 

m(z, z,x) = Y^ M a{z, z)5(x - Q a (z, z)) , 

(33) 

v(z, z,x) = K *m = ^2 M a( z , z)K(x, Q a ) . 

a 



These are peakon solutions for K(x,y) = ^e~\ x ~ y \ . An important feature of these singular 
solutions is that Q a with a £ Z are real functions of z and z. These singular solutions follow 
from Hamilton's principle 5S = for the constrained action S = j L{v : 7, Q) dt given by 



S = J £{u,u) + ^(^M a {z,z)(d z Q a (z,z)-u(Q a ,z,z))\ dsdt. 



Proof. The proof follows the same lines as the proof in j!4] for the real case, keeping in mind 
that x E R is a real number. 

The solution parameters Q a (z, z) € R and M a (z, z) 6 C with a € Z that specify the singular 
solutions (|15p are determined by the following set of evolutionary PDEs in z and z, in which we 



10 



D.D. Holm and R.I. Ivanov 



denote K ab := K(Q a , Q b ) with integer summation indices b, c, e € Z: 
8 z Q a (z, z) = v(Q a , z,z) = J2 M b (z, z)K ab , 

b 

dK ac 

d z M a (z, z) = - d- z M a - Y^(M a M c + M a M c )-^ (no sum on a), (34) 

BK ec 

8 z M a (z, z) = 8- z M a + ^2(M b M c - M b M c ) — (K eb - K cb )(K~ l ) a e ■ 

b,c,e 

The last pair of equations in (|34|) may be solved as a system for the momenta, or Lagrange 
multipliers (M a ), then used in the previous pair to update the support set of positions Q a (z, z). 
Given Q a (z,z) for a € Z, one constructs m(x,z,z) along the solution paths x = Q a (z,z) from 
the first row of (|15p and then obtains u(x,z,z) for x € K from the second row. Alternatively, 
knowing the position Q a (z,z), a € Z, for all z and z, also determines M a upon inverting the 
matrix K ab in the third equation in (|34|) . Note that the third equation in (|34|) is not the complex 
conjugate of the second one. 

Single-peakon solution. The single-peakon solution is again straightforward to obtain from 
equations (JM|) . However, instead of the Laplace equation (TP7|) the solution Q 1 (z,z) satisfies the 
two-dimensional wave equation. 

Two-peakon solution. The two-peakon collision for the complexified G-strand equation 
emerges by following the same path as in the previous section [2^3l with appropriate minor changes. 
In particular, equations (|24p and (|25p become 

Ml = 2(K - K) + 2^' 

M 2 = - ° z - +-^. 
2(K - A") 2 

Now substituting into the second and third equations of the system (|34p and introducing the 
notation 

d z = + iS,, and d s = — id v 
yields the PDE system, 

2(i^o + K) 



^« - * w + 2{K f_ K) ( x t ~ X r,) = ~ \k'{K q - K)(cf>l - <%) . 



(36) 



When <j) = 0, we again have peakon-antipeakon collisions with a simple solution as in section 
12.31 except it involves the two-dimensional wave operator instead of the Laplacian: 

{dl-d 2 v )F(X) = 0, (37) 

where F is the same function as in (|28p . The obvious difference is that the Laplacian d 2 + d% 
in ([25]) for real parameters (s, t) is replaced in ([37]) by the wave operator <9| — d 2 = d z &2 for 
complex parameters (z,z). This type of difference has been studied before in connection with 
the relation between chiral models and harmonic maps in [27] . See [9] for further discussion and 
a guide to the modern literature. The corresponding studies for G-strands can be expected to 
be developed further in future work. For example, it may be interesting to study the extension 
of G-strand equations to supersymmetric variables following methods considered in [10^ 124] , 
particularly for G = Diff(R), if the singular solutions would be preserved. 
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4 Conclusions 

This paper has provided a framework for studying the G-strand equations, comprising the system 
of partial differential equations obtained from the Euler-Poincare (EP) variational equations © 
for a G-invariant Lagrangian, coupled to an auxiliary zero- curvature equation. The latter has 
often been the departure point and main focus in other approaches, especially for the integrable 
chiral models, where it sets up the Lax-pair formulation of the system |18j . 

Once the G-invariant Lagrangian has been specified, the combined dynamics and zero- 
curvature system of G-strand equations in © follows automatically in the EP framework. 

The single-peakon solution of the Diff (R)-strand equations © has been obtained as a har- 
monic function. The stability of the single-peakon solution under perturbations into the full 
solution space of equations © is an interesting problem for future work. 

The dynamics of the singular peakon solutions for the G-strand equations with G = Diff (R) 
has also been studied. The peakon-peakon and peakon-antipeakon collision interactions have 
been found to admit elementary solution methods which again come down to solving Laplace's 
equation for a harmonic function. 

Finally, the complexified version of the Diff G stand equations has been introduced and its 
peakon collision solutions have been formulated and again solved by elementary means. In the 
latter case, the analogous peakon-peakon and peakon-antipeakon collision interactions involved 
the solution of the wave operator, rather than the Laplacian, which occurred in the former case. 

These results lift the previously known relationships between chiral models and harmonic 
maps |27| on finite groups, particularly 0(3), into the EP framework on G = Diff(R). The 
singular solutions (peakons) for both chiral models and harmonic maps on G = Diff (]R) are then 
obtained as momentum maps. 
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